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Abstract 

We divide the cross section for a meson-meson reaction into three parts. The 
first part is for the quark-interchange process, the second for quark-antiquark anni- 
hilation processes and the third for resonant processes. Master rate equations are 
estabhshed to yield time dependence of fugacities of pions, rhos, kaons and vetor 
kaons. The equations include cross sections for inelastic scattering of pions, rhos, 
kaons and vector kaons. Cross sections for quark-interchange-induced reactions, 
that were obtained in a potential model, are parametrized for convenient use. The 
number densities of vr and p (K and K*) are altered by quark-interchange processes 
in equal magnitudes but opposite signs. The master rate equations combined with 
the hydrodynamic equations for longitudinal and transverse expansion are solved 
with many sets of initial meson fugacities. Quark-interchange processes are shown 
to be important in the contribution of the inelastic meson-meson scattering to evo- 
lution of mesonic matter. 
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1 Introduction 



Deconfined matter with high temperature and high density is the focus of studies in 
ultrarelativistic heavy-ion colhsions and its confirmation is well known to be related to 
final-state observables. But hadronic observables are affected by hadronic matter that 
succeeds deconfined matter and measurements on dileptons and photons suffer from a 
background that comes from hadronic matter. In order to clearly identify deconfined 
matter from hadronic observables and electromagnetic probes, one has to subtract any 
influence of hadronic matter. This forces us to pursue a precise description of hadronic 
matter. In addition, a complete knowledge of ultrarelativistic heavy-ion collisions also 
requires understanding of the evolution of hadronic matter. 

Transport models (THIl] can provide with vivid and valid descriptions for the evo- 
lution of hadronic matter. The models deal with known and unknown cross sections 
for hadron-hadron reactions in the following ways. If experimental data are available, 
parametrizations fitted to the data are first used. If no measured data exist, cross sec- 
tions can stem from theoretical calculations, are simply assumed to be constants, are 
adapted to the Breit-Wigner formula with hypothetical widths, or are obtained from the 
parametrizations via the detailed balance or via approximate isospin relations. Features of 
a sort of reactions can be definitely shown by the parametrizations that depend on thresh- 
old energies and the center-of-mass energy of the two colliding hadrons. The assumptions 
of constant cross sections and the approximate isospin relations bring uncertainties to 
study of the evolution of hadronic matter and predictions on final-state observables. 

Time evolution of meson density, for example, meson density [15j, can be studied 
with a rate equation that includes cross sections for scattering of the meson by other 
hadrons. The number of necessary cross sections is huge if rate equations for many 
hadrons are established. To avoid complexity caused by a great number of mesonic degrees 
of freedom in hadronic matter, effective numbers of pions and kaons were introduced to 
approximately get evolution of pionic matter and kaonic matter [TGHTS] . The effective 
number of pions (kaons) is defined as the number of pions (kaons) plus the sum of other 
hadrons weighted by their effective pionic (kaonic) content. For instance, a p meson 
counts as two pions and an Vl baryon counts as three kaons. With this consideration rate 
equations for the effective numbers of pions and kaons were established [T6HT8]. 

At mid-rapidity the PHENIX Collaboration [19] obtained ratios of pT-integrated me- 
son yields in central Au+Au collisions at y^S7Viv=200 GeV as follows: 7r^/7r+ = 0.984, 
K-/K+ = 0.933, K+/TT+ = 0.171 and R-fn- = 0.162, which agree with the results 
of the other collaborations [20] - [22] . The STAR Collaboration found p'^/vr^ = 0.169 in 
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peripheral Au+Au collisions at y/s^jy = 200 GeV [23]. Obviously, pions, rhos and kaons 
are dominant meson species in hadronic matter. In the present work we consider only 
pions, rhos, kaons and vector kaons that constitute mesonic matter. Inelastic meson- 
meson scattering originates from resonances, quark-antiquark annihilation processes and 
quark-interchange processes. At the lowest order a resonance is made up of the remain- 
ing quark and antiquark of which one absorbs a gluon from the annihilation of a quark 
in an initial meson and an antiquark in another initial meson. Of course resonances 
may be glueballs, multiquark states or hybrid states [2ll[25]. At the lowest order the 
quark-antiquark annihilation means that a quark and an antiquark each in a final meson 
come from the annihilation of a quark and an antiquark each in an initial meson. A 
quark-interchange process allows such a type of meson-meson scattering where one gluon 
exchange and the interchange of two quarks each from an initial meson happen. One 
is acquainted with the inelastic meson-meson scattering due to the resonance and the 
quark-antiquark annihilation, but not the scattering induced by the quark interchange. 
One also doesn't know how important the scattering induced by the quark interchange is 
in the evolution of mesonic matter. Indeed, nobody had calculated cross sections for the 
seven quark-interchange-induced reactions tttt ^ pp for 1 = 2, KK -H- K*K* for / = 1, 
KK* o K*K* for I = 1, ttK ^ pK* for / = |, ttK* ^ pK* for J = |, ^ pK* 
for / = I and ttK* ^ pK for / = |, until we obtained cross sections for these reactions 
in Ref. [26] in the quark- interchange mechanism [271|28]. Therefore, we curiously study 
the role of quark-interchange processes that lead to the seven reactions and other isospin 
channels in the evolution of mesonic matter in the present work. The study resorts to 
master rate equations for mesons where reactions of pions, rhos, kaons and vector kaons 
are taken into account. It will be shown that quark- interchange processes are important 
in the contribution of the inelastic meson-meson scattering to the evolution of mesonic 
matter. Therefore, if we include the resonances and the quark-antiquark annihilation 
processes, the quark-interchange processes should be included on an equal footing. 

In the next section the master rate equations for pions, rhos, kaons and vector kaons are 
established while the inelastic meson-meson scattering due to the resonance, the quark- 
antiquark annihilation and the quark interchange is considered. In Section 3 parametriza- 
tions of cross sections for the quark-interchange-induced reactions are presented. Cross 
sections for the quark-antiquark annihilation processes and the resonant processes are 
individually introduced. Numerical results of the master rate equations associated with 
longitudinal expansion and discussions are given in Section 4. The master rate equations 
are extended to include 2^1 mesonic processes in Section 5 and in the case of both the 
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longitudinal and transverse expansion the importance of the quark-interchange processes 
is examined in Section 6. Summary is in the last section. 



2 Master rate equations 

We estabhsh the notation ~ 1 j ^ ~ \ K" pseudoscalar 

isospin doublets as well as K* — I ^ ^ I and K* = \ ^ I for the vector isospin 

doublets. We only consider vr, p, K, K, K* and K* of mcsonic matter. In order to 
clearly exhibit the role of the quark-interchange processes, we neglect both 2 — )■ 1 mesonic 
reactions and decays of p, K* and K* in this section and the next two sections. The 
negligence does not affect us to draw correct conclusions. The inclusion of the 2^1 
mesonic processes is deferred to Section 5. The reactions that can change the numbers of 
7r, p, K, K, K* and K* in a lifetime of mesonic matter are the following inelastic 2-to-2 
scattering: 

1. TTTT -H- pp, 

2. KK o K*K* and KK o k*K*, 

3. KK* ^ K*K* and KK* ^ K*K*, 

4. KK^K*K*, 

5. KK* ^ K*K* and K*K ^ K*K*, 

6. ttK* -h- pK and nK* -H- pK, 

7. nK -H- pK* and nK -H- pK*, 

8. ttK* ^ pK* and tiK* ^ pK*, 

9. pK o pK* and pK o pX*, 

10. TTTT ^ /Ti?, 

11. 7rp O XX* and 7rp O K*K, 

12. KK ^ pp. 
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The cross sections for these reactions are not independent of each other, e.g., cfkk-^k*k* — 

^■kK-^pK* — CTttK-^pK*, CrnK*->pK* — (^■nK*^pK*, CTpK^pK* — CTpK-^pK*, ^nd Cr^p-^K*K — 
'^np-^KK" ■ 

Meson number densities change with time according to the following rate equations, 

d^imu'^) = (1) 

where u'^ = {vP,u) = 7(1, v) is the four-velocity of the local reference frame comoving at 
velocity v and with the Lorentz factor 7. 11,^, n^, n^, Uj^, Uk* and n^, arc the number 
densities of n, p, K, K, K* and K* if i denotes vr, p, K, K, K* and K*, respectively. 
Zero values of the source terms mean that the total number of each particle species is 
conserved. The source terms are given by 



^TT —2 X —{o'pp^TrTrVpp)n^ 2 X — (cr^Tr^ppfTTTr)^^ 

+ {(^pK-^-KK*VpK)npnK — {(y-KK*-^pKVT:K*)n-^nK'' 

+ {'^pK^-kK*'^pk)'<^p1^K ~ {(^ttK*^pK'^ttK*)'I^'k'>^K* 

+ {o'pK*^TTKVpK'-)n'pnK* — {'^TTK^pK*VT,K)n'T,nK 

+ {(^pK'-^-KKVpK*)^p'^K* ~ {'^■KK^pK*V-nK)^'K'^K 

+ {(7pK*-¥'KK*VpK*)npnK* — {(7-KK*-^pK*V-KK*)nT^nK* 

+ {(^ pR'-^-kK*'^ pK*)npnx* — (cr^ft-«_^pK*'*^7r^*)^7r?T'^* 

+ 2{aKK^^^VKK)nKnK - 2 x ^{a^^^KKVnn)nl 

+ {(^KK*-^-Kp'^KK*)'>^K'n'K* — {o'T^p-^KK*'t^Trp)'n'Trnp 

+ {(^K*K^np'"K*K)'nK*nK - {(^np-^K*KV7rp)nnnp, (2) 
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+ {cr7TK*-^pKVnK*)'n.,:nK* — {(y pK^-rrK'VpKjnpnK 

+ {(yiTK-^pK*ViTK)n.,:nK — {(^ pK*-^'KKVpK*)npnK* 
+ {^TrK^pK*'"TrK)'>T'TT'>T'K ^ {^^ pK* ^ttRV pR*)"^ p'^K* 
+ {0'nK*^pK*VnK*)'nnnK* — {o'pK*^TzK*VpK*)npnK* 
+ {'yiTK*->-pK*'t^iTK*)'>^Tvn'K* ~ {(^ pK*-^'KK*'^ pK*)'>^p'>^K* 
+ {'yKK*^iTp'^KK*)'>T'Kn'K* — {o'iTp-^KK*'t^TTp)'>T'Tt'>T'p 
+ {(TK''K^-Kp'^K*K)'>^K*nK — {o'i:p^K*K't^TTp)n'-n-np 

+ '^{crKK^ppVKK)nKnK - 2 x ^{app^KKVpp)nl, 
=2 X ^(cri^*i<r*-s.i<:;<r^;K*K*)^|'* - 2 x ^{aKK^K*K*VKK)nK 

+ -^{<^K*K*-^KK*'VK*K*)n\* — {crKK*-*K*K*VKK*)'nK'nK* 
+ {(TK*K*^KK'^K*K*)'I^K*'nK* — {(TKK^K*K*VKK)'<^KnK 

+ {(J-„K*^pKVT,K')nT,nK* - {<JpK^nK*VpK)npnK 

+ {crpK*-^TrKVpK*)npnK* — {crTrK-^pK*VTrK)nTrriK 

+ {crpK*^pKVpK*)npnK* — {crpK--,pK*VpK)npnK 



+ {'^np-^KK*'t^np)'>T'Trnp — {crKK*-^np'f^KK*)lT'KnK* 
1 

2* 



+ \ {(ypp^KKVpp)nl - {(^KK^ppVKK)nKnK, 
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^K* =2 X -{(JKK->-K*K*VKK)'nK — 2 X -{crK*K*^KKVK*K*)nK* 
+ {cKK*-^K*K*VKK*)'nK'nK* — ■^{cK*K*-^KK*VK*K'')'nl^* 
+ {(^KK^K*K*'^Kk)'<^K'<^K ~ {(^K*K*^KK'^K*K*)'I^K*nK* 

+ {<y pK^nK*VpK)npnK — {o'TTK'^pKVT:K'')nT:nK' 

+ {(y-KK-^pK*VT^K)n-^nK — {crpK*-^-KKVpK*)'npnK* 
+ {crpK-^pK*VpK)npnK — {crpK*-^pKVpK*)npnK* 

The source term of K (K*) is not shown since {"^k*) is obtained from (^E'x*) by 
the replacements of the subscripts, K ^ K and K* ^ K* . Those terms what contain 
the factor 2 relate to the reactions that have two indistinguishable initial or final mesons. 
This factor means that the two initial mesons of a species vanish to attain no final mesons 
of the species or the two final mesons of a species appear from no initial mesons of the 
species. The first fourteen terms of ^^tt equal the negative of the first fourteen terms 
of ^ p. The first fourteen terms of k equal the negative of the first fourteen terms of 
^K* ■ The quark-interchange processes are only contained in these terms. Therefore, the 
quark-interchange processes contribute to the variations of the number densities of tt and 
p, or K and K*, in equal magnitudes but opposite signs. 

The thermal averaged cross section with the relative velocity of two initial mesons v^ei 
is defined as 

where fi{ki) and fj{k2) are the momentum distributions of the two initial mesons with 
the four-momenta ki and k2, respectively; aij^i'j'{y/s) is a cross section that depends 
on the center-of-mass energy ^/s of the two initial mesons. We take the approximate 
factorization form of the Jiittner distribution with nonequilibrium fugacity Aj of particle 
species i, 

m = (7) 

where T is temperature. If two initial mesons are indistinguishable, /j and fj possess the 
same fugacity and as seen in Eqs. (2)- (5) a factor of | is in some terms to remove the 
double counting of initial mesons in the thermal average. 
The number density of particle species i is given by 

d^k A. 



with 

where Wj is the mass of particle species i; the spin-isospin degeneracy factor = 3 for tt, 

— _ — * 

9 for p, 2 for K or K, 6 for K* or X*; k' is the particle momentum in the local comoving 
reference frame. The derivative of hi with respect to T is 

^ = 1 (3n, + , (10) 

with 

For symmetric matter = Xk, A^. = A^*, hx — hx, h^* — hx*, h^- — hx- and 
hx*_ — flK*-- 

Inserting Eqs. (8) and (10) into Eqs. (2)-(5), we obtain rate equations for fugacities 
of TT, p, K and K* of symmetric matter in the longitudinal expansion 

K ( ^ n^-\ T 1 



A,r \ h-^ J T T 

1 Q A^n^ 1 Q 

, o/ \ Q^P^P^KriK V Ox - 

+ 2,{apK^'!TK*VpK}U — 'i{crnK'-^pKVnK*}U Ar'TIk* 

, o/ \ Q^P^P^K*nK'' ^1 \ 0\ - 

+ ^\<7pK''^TxKVpK*jU -—^ Z[aTrK->-pK*V7rK)U AkUk 

, o/ \ O^P^P^K-hK* . > Ox - 

+ ^\0'pK*-^TvK*VpK*)U —— Z\0'TrK*^pK*VTrK*)U AK*TiK* 

Aj^Tlq^ 

\^ I 
+ '^{(7kk-^^^Vkk)u^ - 2 X -((7^^-^;^^v^^)M°A^n^ 

+ 2(cr;^i?*^^pW^K*)M ^{<^^p^KK*v^p}u XpUp, (12) 
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Xp \ rip J T T 

1 oK^^l 1 

— 2 X ~ (c^TT^pp't'TrTr)'?^ T ~ 2 X —{p'pp^T^^Vpf^U Xpflp 

, o/ \ Q^T^ri-K^K*riK* \ Ox - 

+ 2{a^K*-^pKVnK*)U — 2{crpK-^7TK*VpK)U XkUk 

ApTlp 

, X oX^fiT^XKriK „, V o\ - 

+ 2\(y-KK^pK*V-KKlU — 2\CrpK*--fKKVpK*)U AK*nK* 

ApUp 

. ^1 V QXT,nT,XK*nK* „, > On - 

^ ■Z\(y'KK*^pK*VT^K'')U 7-^ Z\(y pK*^T:K*VpK*lU AK*nK* 

ApUp 



, O/ \ O ^K'nKXK*nK* \ Ox - 

+ ^\'^KK*^7!-p'^KK*/'^ ■^\<^iTp^KK''1^np/U At^TIt, 

ArJlr 



^p'^p 

X^.fi^ 1 
+ '^{(^kk^ppVkk)u^^-^ - 2 X -((Tpp^;^;^Vpp)M"A^n^, 

ApTip Z 



Xk ( nx- \ T 1 
3+^ - + - 
Xk \ UK J T T 

1 A2 1 

= 2 X -\(Tk*K*-¥KKVk*K*)u ^ _ 2 X -{aKK-^K*K*VKK}U XkUk 

1 A2 «2 

/ \ A"* X* / \ 0\ — 

+ -\(7k*k*-^kk*vk*k*)u — \(7kk*^k*k*vkk*)u XK*nK* 

2 XrUk 

/ \ K* K* I \ 0\ — 

+ \^k*k''^kk'^k*k*)'^ ~\ — = \Pkk-^k*k*'^kk)'^ XkUk 

XkUk 

/ \ K"* i^* / \ Ox — 

+ \'^K*K*^KK*Vk*K*}'U \'^KK*^K*K*VkK*)U Xk'TIk* 

+ {Cr%K*^pKVT^K*)U r — \(^ pK-^-i:K*V pRjU ApUp 

ARflR 

. I \ oXpnpXK*nK* I \ Ox - 

+ \PpK*--,-KK'^pK*)U 7 — \(7-kK-^pK*V-,^K)U A^riTr 

AkUk 

, / \ 0^pripXK*nK* , \ Ox - 

+ {apK*^pKVpK*)U r-^ {(JpK^pK*VpK)U XpUp 

AkUk 

1 A^n^ 

. . g A,r'^7r Ap'^p ; \ Ox - 

+ Wnp^KK'V^pjU {crKK*^np'^KK*)U XK*nK* 

1 A^n^ 

+ - {cTpp-^KK'Vpp) - ^ - {(^KK-^ppVKK)u^>^KnK, 
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+ 3 + ^^ - + - 



1 1 

= 2 X -{(JkK^K*K*Vkk)u^ . ^ 2 X -{(TK*K*^KKVK*K*)u^XK*nK* 

/ \ K K I \ \ — 

^XpflpXKflK I \ 0\ - 

+ \<^pK^TvK''VpK)U — \0'ttK*^pKVttK*)U ^71^-^ 



+ \<^-kK^pK*V-kK)U — r \0'pK*->-TTKVpK*)U Ap,op 



. I \ QXpUpXxnK I \ 0\ - 

+ KPpK^pK'VpKlU — {crpK*^pKVpK*)u XpUp 

+ {(J^p^K'KV^p)u^ ^^n pnp _ (^aK,K^^^VK*K)u°XKnK, (15) 

where the overdots denote the derivative with respect to the proper time r. 

For the purpose of studying the role of quark-interchange processes, it is enough to only 
consider the longitudinal expansion of hadronic matter. The relativistic hydrodynamic 
equation is 

di^T^" = 0, (16) 
where T'^'' is the energy-momentum tensor given by 

T"'' = {e + P)u^u'' - Pg'"', (17) 

where e is energy density and P is pressure. The simple form of T'^'^ above holds for an 
ideal fluid where viscosity effects have been neglected. The Bjorken's scaling solution of 
the hydrodynamic equation is 

^+i±^ = 0. (18) 
ar T 

The energy density is 

€ = + ep + ex + GR + (^K* + ^R* (19) 

where e^,, e^, ex, €r, ex* and eR* are the energy densities of vr, p, K, K, K* and K* , 
respectively. In order to solve Eq. f lTSj) . relation of pressure and energy density is needed. 
Detailed studies have shown that the relation can be P = 0.15e. Once cross 

sections for the reactions concerned are known, Eqs. (fT2l) - (fT5l) and (fTSjl determine time 
dependence of T(r), Xt,{t), '^pij)i '^k{t) and Xk*{t). 
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3 Cross sections for meson-meson reactions 



The meson-meson cross sections entailed in the master rate equations in Section 2 
are the isospin-averaged cross sections that are obtained by taking the average over the 
isospin states of the two initial mesons and the sum over the isospin states of the two final 
mesons 

a,,^,,,{^s) = (2/^ ^ ^1(2/2 + 1) ^^^^ + 

where Ii and I2 are the isospins of the two initial mesons, respectively, and a{I, y/s) is 
the spin-averaged cross section for the reaction with the total isospin /. 

We explicitly decompose the cross section a {I, y/s) into three parts: the first is ciqi 
from the quark-interchange process, the second (Janni from the annihilation processes and 
the third cTres from the resonant processes. Since the momenta and the coordinates of 
the quark and antiquark constituents of the final mesons in the quark-interchange process 
^{QiQi)+B{q2q2) — > C{qiq2)+D{q2qi) are different from those in the annihilation processes 
or the resonant processes A{qiqi) + 5(^2^2) C{qiqi) + ^(52^2)5 there is no interference 
between the quark-interchange process and the annihilation processes and between the 
quark-interchange process and the resonant processes. No interference of the annihilation 
processes and the resonant processes is usually assumed. Then the cross section for a 
reaction is written as 

y/s) = Cqi(Jqi(/, a/s) + CanniO-anni(-?", V^) + Cj-esC^j-esil , y/s) . (21) 

The coefficients Cqi, Canni and Cres that take values of or 1 are listed in Table [H In the 
following three subsections a^i, cr anni and cxres are presented. For two isospin channels 
all the three processes contribute. For most isospin channels only the quark-interchange 
processes or the annihilation processes contribute. Quark-interchange-induced reactions 
refer to the channels where only the quark interchange works. The quark-interchange 
processes include the quark-interchange-induced reactions. 

3.1 Cross sections for quark- interchange processes 

In Ref. p6] we have obtained unpolarized cross sections for some meson-meson nonreso- 
nant reactions governed only by the quark-interchange mechanism. The cross sections rely 
on mesonic quark- antiquark wave functions and constituent-constituent interaction. The 
quark- antiquark relative-motion wave functions are determined by the Buchmiiller-Tye 
potential [33] that arises from color confinement and one-gluon exchange plus one- and 
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two-loop corrections. The constituent-constituent interaction includes the Buchmiiller- 
Tye potential what is nonrelativistic, central and spin- independent, and spin-spin terms 
that are obtained by performing Foldy-Wouthuysen canonical transformations to a rela- 
tivistic two-constituent Hamiltonian that includes the linear confinement and a relativistic 
one-gluon-exchange potential plus perturbative one- and two-loop corrections [35]. The 
wave functions and the interaction can reproduce the experimental mass splittings be- 
tween the ground-state pseudoscalar octet mesons and the ground-state vector nonet 
mesons [35]. The wave functions and the interaction were employed [26] to calculate cross 
sections for nonresonant reactions A{qiqi) + 5(52^2) — ^ ^(91 ^2) + -D(g2^i) which are en- 
dothermic or exothermic. The A/i-dependence of numerical cross sections exhibited in 
Ref. [26] show a peak for each of endothermic reactions and large magnitudes very near 
the threshold energies of exothermic reactions. For convenient use of the numerical cross 
sections in solving the master rate equations, parametrizations similar to Ref. |36j read 



(Tqi(/, y/s) = CTmax ( — ^ ) exp a ( 
V ^cndo / L V 



(22) 

f cndo , 

where e = ^/s — ^/sq shows difference from the threshold energy Values for the 

parameters cxmax, Cendo and a for the quark-interchange-induced endothermic reactions 
given in the introduction are shown in Table [2l 

The right-hand side of Eq. (22) indicates the zero value of cross section at the threshold 
energy and a fall of cross section at ^/s — )■ 00. Ccndo is close to the center-of-mass energy at 
which an endothermic reaction reaches maximum cross section (Jmax- The power function 
and the exponential function do leave a curve that is not symmetric with respect to the 
peak. The peak is in the energy region that is accessible to meson-meson reactions in 
mesonic matter. 

Table 2 shows the channels with the highest isospins. Since quark-interchange pro- 
cesses can also take place in low isospin channels as seen in Table 1, flavor matrix elements 
/flavor (-^) of the quark-interchange processes are listed in Table |3]for different isospins. The 
entry /flavor (1) = for tttt — t- pp for 1 = 1 means that the reaction in / = 1 is forbidden. 
This is consistent with the fact that the antisymmetric state of vrvr is not allowed. The 
discrepancy of the cross sections for different isospin channels of a reaction results from 
the different flavor matrix elements while matrix elements involving spin and spatial wave 
functions are equal. Then, for instance, we have crqi(/ = 0, ^i) = \<7qi{I = 2, -y/i) and 
(jqi(/ = 1, y/s) = for Tin — pp. Let /max denote the highest isospin of a reaction. The 
cross section for a channel listed in Table 3 is 

/flavor (- ^) 
/flavor (-^maxV 



^qi(A V^) = ,2 Tr \ ^qi(-^ma., V^) • (23) 
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The cross section for an exothermic reaction i'^S^Is) + j^S^I^) — i- i{SiIi) + j{S2l2) is 
obtained from the endothermic reaction ij — )■ i'j' by the detailed balance 

where = (25i + 1) (2/i + l) (2^2 + 1) (2/2 + 1) and = (2^3 + 1) (2/3 + 1) (2^4 + 1) (2/4 + 1) 
denote the spin-isospin degeneracy factors of initial particles with the spins Si and 5*2 as 
well as the isospins Ji and I2, and of final particles with the spins 6*3 and 5*4 as well as 
the isospins J3 and I4, respectively; P and P' denote momenta of an initial meson and a 
final meson in the center-of-momentum frame of the reaction ij — i'j', respectively. 



3.2 Cross sections for annihilation processes 

The isospin-averaged cross section for an annihilation reaction can be parametrized 
by [371138] 

0-anni(v^)=6(l-^)'. (25) 

Values of the parameter b and the dimensionless parameter c for the reactions nn — )■ KK, 
vrp — )■ KK*{K*K) and KK — )■ pp are given in Table HI Since no experimental data are 
available, we assume that cross sections for other annihilation processes listed in Table 
1 possess Eq. (25) with the same parameters h and c as the reaction vrvr — )■ KK . The 
treatment of the annihilation processes is simple. 



3.3 Cross sections for resonant processes 

Cross sections for resonant processes are generally described by the Breit-Wigner for- 
mula 

/ j-s^ _ 2J + 1 vr r^i?ini?out 

''■"'^^^^'^ " (25i + 1)(252 + 1) (v^ - mnP + T^/A ' ^^^^ 

where a resonance has its spin J, its energy and its full width F, and B-^^ and i?out are 
the branching fractions of the resonance decays into the initial state and the final state, 
respectively. We take account of the resonances /o(1370), p(1450), /o(1500) and p3(1690) 
for TTTT o pp, fs:i(1270), Ki{Um), is:* (1410), 7^2* (1430), fs:*(1680) and is:3*(1780) for 
TiK* ^ pK, /o(980), 0(1020), /2(1270), /o(1370), p(1450), /o(1500), /^(1525), P3(1690), 
p(1700), /o(1710), /2(1810) and ^(2050) for nn ^ KK. Since the full widths and the 
branching fractions for some resonances have not been fixed by measurements [321, their 
values we select for the reactions tttt •(-)■ pp, ttK* -h- pK and tttt -h- KK are listed in Table 

El 
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4 Results and discussions 



In this section we represent and discuss results that are from solving the master rate 
equations in combination with the hydrodynamic equation (18) simultaneously by numer- 
ical integration using a fourth order Runge-Kutta method. The inelastic 2-to-2 scattering 
in the master rate equations includes the three types of processes: the quark-interchange 
processes, the annihilation processes and the resonant processes. As a good approxi- 
mation, we assume that hadronization of quark-gluon plasma at the critical temperature 
Tc = 175 MeV [10] only produces vr, p, K, K* and K*. We assume that the hadroniza- 
tion is finished at = 5.6 fm/c and at the moment mesonic matter has A,r = 0.7, \p = 0.7, 
Xk = 0.5 and A^-* = 0.2 for the fugacities of vr, p, K and K*, respectively. We consider 
mesonic matter at or near mid-rapidity, i.e., m'^ ~ (1, 0, 0, 0), solve Eqs. (12)-(15) and (18), 
and terminate numerical calculations when mesonic matter reaches the kinetic freeze-out 
temperature Tf^ = 105 MeV, which corresponds to a freeze-out time of the order of 30 
fm/c. 

If the inelastic 2-to-2 scattering is switched off, ie., the source terms are zero, the 
master rate equations become 



Together with Eq. (8) we have 



For massive bosons, nj given by Eq. (9) and energy densities do not simply rely on a power 
of T, and the hydrodynamic equation cannot guarantee UiT as constants. Therefore, Aj 
depend on the proper time r unlike the case of massless bosons which fugacities are 
constants. Fugacities for vr, p, K and K* are denoted by X-^no, Apno, Axno and XK*no, 
respectively, and are plotted as solid curves in Fig. 1. The solid curves for p, K and K* 
rise with the increase of time. The fugacity of the lightest meson first decreases slightly 
and then increases. Compared to the results in the absence of the source terms, we 
show meson fugacities by the dashed curves and indicate the meson fugacities by Ajqar 
[i = 7c, p, K, K*) while the quark-interchange processes, the annihilation processes and 
the resonant processes are all included. The differences between Ajqar and Ajno due to 
the inelastic 2-to-2 scattering are obvious. To show how the quark-interchange processes 
modify fugacities, we show the fugacities by dotted curves and denote the fugacities by 
Aiqi {i = 71, p, K, K*) while only the quark-interchange processes govern time dependence 




(27) 



which have the solutions 
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of the fugacities. In most of the range 5.6 fm/c < r < 30 fm/c the absolute values of the 
fugacity differences, | Ajqi — Ajno | , caused by the quark-interchange processes are smaller 
than I Ajqar — Ajno | causcd by the three types of processes. 

To show a role of the quark-interchange processes, we define 

Ri = ^'q' ~ . (30) 

The larger the absolute values of the more important the quark-interchange processes. 
If Ri > 0, either of the quark-interchange processes and the combination of the three types 
of processes increases (reduces) fugacities relative to Ajno- If Ri < 0, the quark-interchange 
processes increase (reduce) fugacities relative to Ajno while the annihilation and resonant 
processes reduce (increase) fugacities. Values of Ri change with the increase of time. At 
T = 20 fm/c, R^ = 0.37, Rp = 0.30, Rk = 0.21 and Rk* = 0.27. 

To quantitatively determine the importance of the quark-interchange processes, we 
define the average of the absolute value of Ri by 

Ri = ^-^ -■ (31) 

where Tfz is the freeze-out time of mesonic matter. Then, R^^ = 0.53, Rp = 0.30, Rk = 0.21 
and Rk* = 0.27. Hence, with the set of initial fugacities, A^r — 0.7, Xp = 0.7, Xk = 0.5 
and Xk* ~ 0.2, for the master rate equations and the hydrodynamic equation, the quark- 
interchange processes are important in the contribution of the inelastic 2-to-2 scattering 
to the evolution of mesonic matter. But the conclusion is only limited to this set. We 
need to examine Ri versus other initial fugacities. Initial fugacities of mesonic matter 
depend on incident energies in nucleus-nucleus coUisions. The Au-Au collisions have been 
carried out at various energies of per pair of colliding nucleons allowed by the Relativistic 
Heavy Ion Collider (RHIC) and Pb-Pb coUisions at the Large Hadron CoUider (LHC) 
have been performed at higher energies. Different nucleus-nucleus collisions at different 
energies produce mesonic matter with different magnitudes of initial fugacities. Therefore, 
we use a wide range of initial fugacities to check the importance of the quark-interchange 
processes. It is impossible to plot graphs for Ri {i — tt, p, K, K*) versus the four variables 
Att, Ap, Xk and A^*, but we can tabulate Ri {i — tt, p, K, K*) versus a finite sets of initial 
fugacities in the range between and 1. We take 81 sets in which A^ = 0.35, 0.65, 0.95, 
Xp = 0.15, 0.45, 0.75, A^ = 0.25, 0.55, 0.85 and A^* = 0.15, 0.55, 0.95. Statistically, the 81 
sets can tell us how important the quark-interchange processes are. The freeze-out time 
Tfz depends on initial fugacities. Averages Ri {i = tt, p, K, K*) for the 81 sets of initial 
fugacities are listed in the middle four columns in Tables 6-8. 
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The average may be as large as 33.07 and this occurs when the contributions of the 
quark-interchange, annihilation and resonant processes can cancel each other. The average 
may be as small as 0.01. Most of the entries in the middle four columns are above 
0.2. Therefore, the quark-interchange processes are important in the contribution of the 
inelastic 2-to-2 scattering to the evolution of mesonic matter. 

Via HIJING Monte Carlo simulations [lTI - 115] . initial fugacities of deconfined gluons 
and quarks produced in central Au-Au collisions at a/sJvw" = 200 GeV are about 0.2 and 
0.032, respectively, as seen in the second set of initial conditions of quark-gluon plasma in 
Table I of Ref. [H]. Using the ratio 0.2 of strange-quark number to up-quark number ^5], 
that may reproduce the measured ratios j-n^ = 0.15±0.02 and K*^/K~ = 0.205±0.033 
at midrapidity |16] , we solve master rate equations of quark-gluon plasma given in Ref. [17] 
to obtain time dependence of fugacities of gluons and quarks. We obtain the time Th ~ 5.6 
fm/c, the fugacities Xg ~ 0.75 and Xg ^ 0.52 at T^. Coalescence of quarks and antiquarks 
forms mesons at T^. Assume that the formed mesons are only vr, p, K, K, K* and K* and 
that = Xp and Xk = Xk* at Tc. Then we obtain = Ap ~ 1.31 and Xk = Xk* ~ 0.56 
that are the initial fugacities of mesonic matter produced in central Au-Au collisions at 
^ysNN = 200 GeV. For other nucleus-nucleus collisions meson fugacities are different from 
these values. The pion fugacity may not equal the rho fugacity at and the kaon fugacity 
may not equal the vector kaon fugacity. If ^/sNN decreases continuously, meson fugacities 
decrease continuously. Therefore, A^ = Ap = 0.7, Xk = 0.5 and Aa'* = 0.2 used in this 
section and some sets of initial fugacities listed in Tables 6-8 can be covered in some 
nucleus-nucleus collisions at some values of ^/snn- Since the strange-quark number is less 
than half the up-quark number, those sets of initial fugacities in which both Xk and A^-. 
are larger than A^^ and Ap are not possible. But such sets of initial fugacities yield Ri > 0.2 
{i = IT, p, K, K*) in more than half the entries of Ri. Therefore, the impossible sets of 
initial fugacities help us more firmly establish that the quark-interchange processes are 
important in the contribution of the inelastic 2-to-2 scattering to the evolution of mesonic 
matter. 

It is shown in Table 1 that fourteen reaction channels involve the quark-interchange 
processes, sixteen channels the annihilation processes and three channels the resonant 
processes. In a reaction where a quark-interchange process occurs, the channel with the 
highest isospin is only induced by the quark-interchange process and the number 2/ -|- 1 
of isospin component in Eq. (20) can enhance the capability of the process in influencing 
the evolution of mesonic matter. 
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5 Results pertinent to 2 ^ 1 processes 



Since we only consider vr, p, K, K, K* and K* , resonances involved are p, K* and K* . 
Then we include vrvr -H- p, vri^ ^ i^* and vr.ft' ^ K* . Let rp_j.^7r, Fii-.^Trj^ and F/f*^^^ 
be the decay widths of p — )■ tttt, fC* — )■ Tri^ and K* — )■ vr^, respectively. We add the 
following four expressions 

and 

— ri<-*_j.7r_fs:ni^. + {0'TTK^K*VnK)nTrnK, 

to the source terms \l/p, \E'i^ and \E'it. in Eqs. (2)-(5), respectively, to establish 
master rate equations with the 2 o 1 mesonic processes. Obtained from the experimental 
data [48j, the cross section for tttt — p is 

_ 80 mb 

^^^"'^ " i + 4(yi-m,)vr^-...- ^^^^ 

The cross section for ttK — ?► K* or tcK — K* can be found in Refs. [3111]. From solutions 
of the master rate equations with the 2^1 processes and the hydrodynamic equation 
(18), we get the average values Ri at various initial fugacities and list them in the right 
four columns in Tables 6-8. 

In the tables about 96% of the entries in the right four columns have values larger than 
0.5 and 53% larger than 1. Moreover, most of Ri {i = vr, p, K, K*) are larger than the 
corresponding ones derived from the master rate equations without the 2^1 processes. 
Therefore, while the quark-antiquark annihilation processes and the resonant processes 
are taken into account, the quark- interchange processes must be included on an equal 
footing. 



6 Results pertinent to transverse expansion 

In the preceding sections we have considered only the longitudinal expansion for 
mesonic matter. In this section we rely on both longitudinal and transverse expansion 
to deal with mesonic matter produced in central collisions. The four- velocity of the local 
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reference frame is u'^ = 7(^,fr,0, ^) with 7 = 1/a/1 — f^, where Vr is the transverse 
velocity. The left-hand side in Eq. (1) becomes 

Ofii O'y 7 1.0 
df,{niU^) = 7— + ni{— + -) + - — {rnavr). (33) 
or or T r or 

For matter uniformly distributed, 
In terms of fugacities, 

Oi,[n.iU^) = v^na- — Vu^Xa—- — V2niYVr^ — hni7(27 + — \ \ . (35) 

OT oT OT OT or T r 

Solving the master rate equations with Eq. (35) and the 2^1 mesonic processes and the 
hydrodynamic equations describing the longitudinal and transverse expansion in Ref. ^9] . 
we obtain, for example, at = 0.7, Xp = 0.7, Xk = 0.5 and Xk* = 0.2, the average values 
Ri {i = n, p, K, K*) are 1.11, 0.78, 1.21 and 0.85, respectively. Corresponding to most 
sets of the initial fugacities listed in Tables 6-8, R,,, Rp, Rk and Rk* are larger than 
1. Therefore, we must use the quark-interchange processes on an equal footing while the 
annihilation processes and the resonant processes are considered. 



7 Summary 

We have established a set of master rate equations that describe time dependence 
of fugacities of pions, rhos, kaons and vector kaons in mesonic matter. A meson-meson 
reaction is comprised of the quark-interchange process, the annihilation processes and 
the resonant processes. The cross sections for the quark-interchange-induced reactions, 
that were obtained from the Buchmiiller-Tye potential plus the spin-spin interaction, are 
parametrized for convenient use in studying the evolution of mesonic matter. 

The variations of fugacities of pions, rhos, kaons and vector kaons are governed by 
the inelastic meson-meson scattering, the 2 f-)- 1 mesonic processes and the expansion 
of mesonic matter. In most reactions the quark-interchange processes take place. If 
the number density of vr is increased (reduced) by the quark-interchange processes, the 
number density of p is reduced (increased) in the same amount. This relation also holds 
true for K and K* . Numerical results of the master rate equations show that the quark- 
interchange processes are important in the contribution of the inelastic 2-to-2 scattering 
to the evolution of mesonic matter. 



18 



Acknowledgments 

This work was supported by National Natural Science Foundation of China under 
Grant No. 10675079. 



19 



References 



H. Sorge, H. Stocker, W. Greiner, Ann. Phys. 192, 266 (1989). 

H. Sorge, Phys. Rev. C 52, 3291 (1995). 

B.-A. Li, CM. Ko, Phys. Rev. C 52, 2037 (1995). 

Z.-W. Lin, CM. Ko, B.-A. Li, B. Zhang, S. Pal, Phys. Rev. C 72, 064901 (2005). 

W. Gassing, E.L. Bratkovskaya, Phys. Rep. 308, 65 (1999). 

S.A. Bass et ai, Prog. Part. Nucl. Phys. 41, 255 (1998). 

G. Nonaka, S.A. Bass, Phys. Rev. G 75, 014902 (2007). 

D.E. Kahana, S.H. Kahana, Phys. Rev. G 58, 3574 (1998). 

D.E. Kahana, S.H. Kahana, Phys. Rev. G 59, 1651 (1999). 

B.-H. Sa, A. Tai, Comput. Phys. Commun. 90, 121 (1995). 

A. Tai, B.-H. Sa, Gomput. Phys. Gommun. 116, 353 (1999). 
T.J. Humanic, Phys. Rev. G 57, 866 (1998). 

Y. Nara, N. Otuka, A. Ohnishi, K. Niita, S. Ghiba, Phys. Rev. G 61, 024901 (1999). 

T. Hirano, U. Heinz, D. Kharzeev, R. Lacey, Y. Nara, Phys. Lett. B 636, 299 (2006). 

L. Alvarez-Ruso, V. Koch, Phys. Rev. C 65, 054901 (2002). 

S. Gavin, NucL Phys. B 351, 561 (1991). 

G. Song, V. Koch, Phys. Rev. G 55, 3026 (1997). 

S. Pratt, K. Haghn, Phys. Rev. G 59, 3304 (1999). 

S.S. Adler et al. (PHENIX GoUaboration) , Phys. Rev. G 69, 034909 (2004). 

G. Adler et al. (STAR Collaboration), Phys. Lett. B 595, 143 (2004). 

I. Arsene et al. (BRAHMS Collaboration), Nucl. Phys. A 757, 1 (2005). 

B. B. Back et al. (PHOBOS Collaboration), Nucl. Phys. A 757, 28 (2005). 
J. Adams et al. (STAR Collaboration), Phys. Rev. Lett. 92, 092301 (2004). 

C. Amsler, N.A. Tornqvist, Phys. Rep. 389, 61 (2004). 

D. V. Bugg, Phys. Rep. 397, 257 (2004). 

Y.-Q. Li, X.-M. Xu, Nucl. Phys. A 794, 210 (2007). 
T. Barnes, E.S. Swanson, Phys. Rev. D 46, 131 (1992). 

E. S. Swanson, Ann. Phys. 220, 73 (1992). 
J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

H. Sorge, Phys. Lett. B 402, 251 (1997). 

20 



L.V. Bravina et al, Phys. Rev. C 60, 024904 (1999). 
P.P. Kolb, J. Sollfrank, U. Heinz, Phys. Lett. B 459, 667 (1999). 
P.P. Kolb, J. Sollfrank, U. Heinz, Phys. Rev. C 62, 054909 (2000). 
W. Buchmliller, S.-H.H. Tye, Phys. Rev. D 24, 132 (1981). 
X.-M. Xu, Nucl. Phys. A 697, 825 (2002). 

T. Barnes, E.S. Swanson, C.-Y. Wong, X.-M. Xu, Phys. Rev. C 68, 014903 (2003). 

G.E. Brown, CM. Ko, Z.G. Wu, L.H. Xia, Phys. Rev. C 43, 1881 (1991). 

W. Gassing, E.L. Bratkovskaya, U. Mosel, S. Teis, A. Sibirtsev, Nucl. Phys. A 614, 
415 (1997). 

Particle Data Group, W.-M. Yao et al, J. Phys. G 33, 1 (2006). 
F. Karsch, E. Laermann, A. Peikert, Nucl. Phys. B 605, 579 (2001). 
X.-N. Wang, M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 
M. Gyulassy, X.-N. Wang, Comput. Phys. Gommun. 83, 307 (1994). 
X.-N. Wang, Phys. Rep. 280, 287 (1997). 

X.-M. Xu, D. Kharzeev, H. Satz, X.-N. Wang, Phys. Rev. G 53, 3051 (1996). 
T.S. Biro, P. Levai, J. Zimanyi, Phys. Rev. G 59, 1574 (1999). 



O. Barannikova, for the STAR Gollaboration, arXiv:nucl-ex/0403014vl. 

P. Levai, B. Mliller, X.-N. Wang, Phys. Rev. G 51, 3326 (1995). 

V. Flaminio, W.G. Moorhead, D.R.O. Morrison, N. Rivoire, GERN, Geneva Report 
No. GERN-HERA-84-01, 1984. 

[49] H. von Gersdorff, L. McLerran, M. Kataja, P.V. Ruuskanen, Phys. Rev. D 34, 794 
(1986). 



21 



Table 1: Values of Cqi, Canni and Cj 



Channel 


Cqi 


Canni 


Cres 


I = 2 TTTT -H- pp 


1 








I = 1 TTTT <^ pp 





1 





I = nn ^ pp 


1 


1 


1 


I = 1 KK ^ K*K* 


1 








1 = KK ^ K*K* 


1 








I = 1 KK* ^ K*K* 


1 








1 = KK* K*K* 


1 








I = 1 KK ^ K*K* 





1 





1 = KK ^ K*K* 





1 





I = 1 KK* ^ K*K* 





1 





1 = KK* ^ K*K* 





1 





I = 3/2 ttK* o pK 


1 








1 = 1/2 ttK* ^pK 


1 


1 


1 


I = 3/2 ttK o pK* 


1 








1 = 1/2 ttK^ pK* 


1 


1 





I = 3/2 ttK* o pK* 


1 








1 = 1/2 ttK* ^pK* 


1 


1 





I = 3/2 pK o pK* 


1 








1 = 1/2 pK^pK* 


1 


1 





7 = 1 TTTT -H- KK 





1 





7 = TTTT O KK 





1 


1 


I = l7rp^ KK* 





1 





7 = 7rp O KK* 





1 





I=1KK ^ pp 





1 





7 = KK ^ pp 





1 
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Table 2: Parameters in Eq. (22). 



Channel 


CTmax (nib) 


Cendo (GeV) 


a 


I — 2 TTTT ^ pp 


0.49991 


0.20909 


0.87446 


I = 1 KK K*K* 


0.61622 


0.16539 


0.4883 


I = 1 KK* K*K* 


0.85168 


0.26399 


1.05175 


I = 3/2 ttK* pK 


1.40233 


0.15023 


1.07478 


I = 3/2 ttK pK* 


0.49839 


0.12056 


0.40939 


I = 3/2 nK* pK* 


0.49 


0.21 


0.88 


I = 3/2 pK pK* 


0.5081 


0.3166 


1.89693 


Table 3: Flavor matrix elements /flavor (-^)- 


7 = 


7 = 1 7 


= 2 7 = 1/2 


7 = 3/2 


TTTT — >■ pp ~l/2 





1 




KK K*K* 1 


1 






KK* K*K* 1 


1 






nK* pK 




-1/2 


1 


ttK pK* 




-1/2 


1 


nK* pK* 




-1/2 


1 


pK pK* 




-1/2 


1 



Table 4: Values of b and c. 



Reaction 


b (mb) 


c 


TTTT — >■ KK 


2.7 


0.76 


np KK*{K*K) 


0.4 


0.5 


KK -> pp 


3.5 


0.38 
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Table 5: Some resonances formed in tttt •H- pp, ttK* •h- pK and tttt •H- KK. 



Name 


r (MeV) 




Bpp 


/o(1370) 


370 


0.26 


0.208 


p(1450) 


147 


0.0672 


0.02 


Name 


r (MeV) 




BpK 


i^i(1400) 


174 


0.94 


0.3 


K*{U10) 


232 


0.4 


0.07 


K*{1680) 


322 


0.299 


0.314 


Name 


r (MeV) 




Brk 


/o(980) 


70 


0.755 


0.245 


/o(1370) 


370 


0.203 


0.35 


p(1450) 


147 


0.0672 


0.0016 


p(1700) 


250 


0.2345 


0.0412 


/2(1810) 


197 


0.0048 


0.003 
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Table 6: R-,^, Rp, Rk and Rk* irrelevant and relevant to the 2 -H- 1 processes are listed 
in the middle and right four columns, respectively. Initial fugacities are in the left four 
columns. 





Xp 


\k 


Aa"* 


Rtt 


Ro 


Rk 


Rk* 


R-K 


Ro 


Rk 


Rk* 


0.35 


0.15 


0.25 


0.15 


0.43 


0.18 


1.14 


1.80 


1.34 


1.54 


0.95 


1.02 


0.35 


0.15 


0.25 


0.55 


1.63 


0.19 


0.30 


0.20 


0.99 


2.03 


4.79 


0.93 


0.35 


0.15 


0.25 


0.95 


0.36 


0.19 


0.27 


0.34 


0.87 


2.81 


1.17 


0.90 


0.35 


0.15 


0.55 


0.15 


1.19 


0.22 


0.01 


2.52 


1.75 


6.89 


0.53 


1.06 


0.35 


0.15 


0.55 


0.55 


0.45 


0.23 


0.17 


0.03 


1.34 


2.68 


0.36 


0.85 


0.35 


0.15 


0.55 


0.95 


0.27 


0.21 


2.86 


0.04 


0.93 


2.04 


6.22 


0.83 


0.35 


0.15 


0.85 


0.15 


5.28 


0.50 


0.01 


2.56 


4.49 


2.42 


0.50 


1.72 


0.35 


0.15 


0.85 


0.55 


0.54 


0.74 


0.04 


0.08 


2.23 


2.71 


0.35 


0.77 


0.35 


0.15 


0.85 


0.95 


0.29 


0.22 


0.09 


0.03 


2.34 


4.07 


0.27 


0.77 


0.35 


0.45 


0.25 


0.15 


0.61 


0.28 


0.22 


0.18 


1.07 


1.00 


0.98 


1.16 


0.35 


0.45 


0.25 


0.55 


0.44 


2.76 


0.27 


0.39 


1.00 


5.32 


1.66 


0.96 


0.35 


0.45 


0.25 


0.95 


0.40 


2.16 


0.22 


0.42 


0.93 


1.41 


1.52 


0.94 


0.35 


0.45 


0.55 


0.15 


0.28 


0.21 


0.13 


2.07 


1.27 


0.97 


0.77 


1.77 


0.35 


0.45 


0.55 


0.55 


0.74 


0.26 


0.10 


0.13 


1.48 


1.04 


0.67 


0.93 


0.35 


0.45 


0.55 


0.95 


0.50 


0.66 


0.08 


0.22 


0.99 


1.29 


0.45 


0.90 


0.35 


0.45 


0.85 


0.15 


0.45 


0.20 


0.10 


5.44 


1.98 


0.97 


0.68 


0.80 


0.35 


0.45 


0.85 


0.55 


0.26 


0.21 


0.08 


0.10 


1.98 


1.01 


0.60 


0.87 


u.oo 


n /I 


U.oO 


u.yo 


n QQ 
u.yy 




U.U ( 


nil 

U. ii 


i. ( 


i . io 


u.ou 


U.oO 


0.35 


0.75 


0.25 


0.15 


7.50 


0.28 


0.25 


0.40 


1.04 


0.99 


1.00 


1.25 


0.35 


0.75 


0.25 


0.55 


0.42 


0.41 


0.34 


0.42 


1.00 


1.02 


1.38 


0.97 


0.35 


0.75 


0.25 


0.95 


0.42 


0.57 


0.58 


0.43 


0.95 


1.07 


1.31 


0.95 


0.35 


0.75 


0.55 


0.15 


0.20 


0.24 


0.18 


0.47 


1.16 


0.96 


0.83 


1.72 


0.35 


0.75 


0.55 


0.55 


0.92 


0.28 


0.19 


0.22 


1.12 


0.98 


0.78 


0.96 


0.35 


0.75 


0.55 


0.95 


1.32 


0.34 


0.19 


0.30 


1.03 


1.01 


0.70 


0.93 


0.35 


0.75 


0.85 


0.15 


0.32 


0.22 


0.15 


0.69 


1.45 


0.94 


0.74 


1.11 


0.35 


0.75 


0.85 


0.55 


0.39 


0.25 


0.15 


0.56 


1.42 


0.96 


0.70 


1.11 


0.35 


0.75 


0.85 


0.95 


4.53 


0.28 


0.16 


0.19 


1.72 


0.98 


0.64 


0.89 
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Table 7: The same as Table 6. 





K 


Xk 


Xk* 




Rp 


Rk 


Rk* 


R-K 


Rp 


Rk 


Rk* 


0.65 


0.15 


0.25 


0.15 


0.15 


0.17 


0.14 


1.74 


0.68 


1.44 


0.70 


1.08 


0.65 


0.15 


0.25 


0.55 


0.17 


0.20 


0.20 


1.61 


1.25 


0.73 


0.89 


0.87 


0.65 


0.15 


0.25 


0.95 


0.11 


0.20 


0.22 


1.36 


0.98 


0.74 


0.94 


0.83 


0.65 


0.15 


0.55 


0.15 


10.48 


0.16 


0.14 


5.73 


0.92 


2.30 


7.31 


1.11 


0.65 


0.15 


0.55 


0.55 


0.51 


0.21 


6.34 


0.07 


0.86 


0.68 


5.43 


1.00 


0.65 


0.15 


0.55 


0.95 


2.90 


0.21 


0.51 


0.08 


0.87 


0.72 


5.80 


0.76 


0.65 


0.15 


0.85 


0.15 


1.50 


0.18 


0.03 


4.62 


1.34 


1.42 


0.34 


1.74 


0.65 


0.15 


0.85 


0.55 


2.15 


0.21 


0.24 


0.11 


0.78 


1.09 


0.41 


1.35 


0.65 


0.15 


0.85 


0.95 


1.71 


0.22 


3.52 


0.09 


0.83 


0.69 


12.40 


0.69 


0.65 


0.45 


0.25 


0.15 


0.14 


2.44 


2.00 


0.09 


1.73 


1.67 


3.07 


1.43 


0.65 


0.45 


0.25 


0.55 


0.59 


0.13 


0.18 


0.69 


1.37 


1.64 


0.74 


0.96 


0.65 


0.45 


0.25 


0.95 


1.03 


0.15 


0.21 


0.62 


0.98 


1.99 


0.86 


0.93 


0.65 


0.45 


0.55 


0.15 


0.36 


0.28 


0.13 


0.33 


1.62 


1.04 


0.72 


0.92 


0.65 


0.45 


0.55 


0.55 


2.09 


1.14 


0.80 


0.13 


1.52 


1.60 


1.08 


0.90 


0.65 


0.45 


0.55 


0.95 


0.44 


0.17 


3.52 


0.22 


0.97 


2.56 


2.85 


0.87 


0.65 


0.45 


0.85 


0.15 


7.58 


0.23 


0.09 


0.51 


1.62 


0.97 


0.60 


1.19 


0.65 


0.45 


0.85 


0.55 


0.37 


4.12 


0.03 


0.47 


4.25 


1.53 


0.48 


1.10 


0.65 


0.45 


0.85 


0.95 


0.29 


0.30 


0.08 


0.07 


1.68 


1.91 


0.35 


0.80 


0.65 


0.75 


0.25 


0.15 


1.85 


0.52 


2.24 


0.12 


1.18 


1.05 


3.80 


42.57 


0.65 


0.75 


0.25 


0.55 


0.76 


3.09 


0.10 


0.60 


1.06 


1.46 


1.99 


0.98 


0.65 


0.75 


0.25 


0.95 


0.54 


4.69 


0.19 


0.53 


0.98 


1.85 


0.89 


0.95 


0.65 


0.75 


0.55 


0.15 


0.27 


0.27 


0.20 


0.30 


1.27 


0.97 


0.84 


1.13 


0.65 


0.75 


0.55 


0.55 


0.93 


0.57 


0.16 


0.46 


1.15 


1.09 


0.78 


1.01 


0.65 


0.75 


0.55 


0.95 


0.44 


1.58 


0.51 


0.33 


0.98 


1.90 


0.60 


0.92 


0.65 


0.75 


0.85 


0.15 


0.37 


0.24 


0.15 


0.42 


1.65 


0.94 


0.71 


0.85 


0.65 


0.75 


0.85 


0.55 


0.26 


0.26 


0.13 


0.82 


1.51 


0.99 


0.65 


1.71 


0.65 


0.75 


0.85 


0.95 


0.88 


0.45 


0.10 


0.17 


1.19 


1.16 


0.54 


0.87 
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Table 8: The same as Table 6. 





K 


Ax 


Xk* 




Rp 


Rk 


Rk* 


R-K 


Rp 


Rk 


Rk* 


0.95 


0.15 


0.25 


0.15 


0.12 


0.17 


0.09 


3.27 


0.86 


0.88 


0.65 


1.15 


0.95 


0.15 


0.25 


0.55 


0.16 


0.20 


0.17 


2.74 


0.85 


0.84 


0.78 


1.44 


0.95 


0.15 


0.25 


0.95 


0.17 


0.20 


0.20 


1.32 


1.26 


0.83 


0.85 


1.30 


0.95 


0.15 


0.55 


0.15 


0.23 


0.16 


3.95 


0.34 


0.92 


0.90 


4.40 


0.96 


0.95 


0.15 


0.55 


0.55 


0.26 


0.21 


0.34 


0.08 


0.95 


0.86 


1.38 


1.64 


0.95 


0.15 


0.55 


0.95 


0.28 


0.21 


0.31 


0.21 


1.75 


0.84 


1.29 


1.88 


0.95 


0.15 


0.85 


0.15 


4.44 


0.14 


0.11 


0.95 


0.95 


0.90 


1.54 


1.03 


0.95 


0.15 


0.85 


0.55 


0.67 


0.21 


3.62 


0.10 


1.02 


0.86 


8.08 


1.67 


0.95 


0.15 


0.85 


0.95 


0.73 


0.22 


8.17 


0.11 


1.98 


0.84 


8.36 


2.43 


0.95 


0.45 


0.25 


0.15 


0.06 


0.08 


0.03 


0.03 


2.55 


1.32 


0.30 


1.26 


0.95 


0.45 


0.25 


0.55 


0.06 


0.17 


0.16 


33.07 


1.73 


1.41 


0.65 


0.94 


0.95 


0.45 


0.25 


0.95 


0.04 


0.18 


0.19 


1.84 


1.31 


2.17 


0.77 


0.86 


0.95 


0.45 


0.55 


0.15 


0.14 


1.42 


0.45 


0.24 


1.74 


1.25 


1.15 


0.89 


0.95 


0.45 


0.55 


0.55 


0.15 


0.19 


0.54 


0.34 


1.13 


1.21 


2.25 


1.03 


0.95 


0.45 


0.55 


0.95 


0.10 


0.21 


0.33 


0.28 


1.05 


1.13 


3.67 


0.81 


0.95 


0.45 


0.85 


0.15 


0.97 


0.75 


0.08 


0.37 


1.15 


1.23 


0.47 


0.96 


0.95 


0.45 


0.85 


0.55 


0.99 


0.31 


0.19 


1.44 


1.02 


1.47 


0.57 


1.22 


0.95 


0.45 


0.85 


0.95 


2.48 


0.24 


2.60 


0.05 


0.82 


1.21 


15.76 


1.01 


0.95 


0.75 


0.25 


0.15 


0.17 


5.86 


0.20 


0.08 


2.13 


2.59 


4.36 


1.46 


0.95 


0.75 


0.25 


0.55 


0.46 


0.08 


0.14 


1.11 


1.76 


2.07 


0.41 


0.99 


0.95 


0.75 


0.25 


0.95 


15.83 


0.12 


0.19 


0.72 


1.49 


2.32 


0.68 


0.93 


0.95 


0.75 


0.55 


0.15 


0.96 


0.41 


0.29 


0.25 


1.90 


1.09 


0.93 


1.12 


0.95 


0.75 


0.55 


0.55 


4.09 


4.11 


0.28 


1.85 


1.51 


1.40 


0.85 


1.16 


0.95 


0.75 


0.55 


0.95 


8.19 


0.14 


1.61 


0.41 


1.13 


1.58 


1.39 


0.88 


0.95 


0.75 


0.85 


0.15 


0.32 


0.28 


0.16 


0.34 


1.57 


0.88 


0.66 


0.88 


0.95 


0.75 


0.85 


0.55 


0.25 


0.78 


0.09 


1.51 


1.31 


1.16 


0.54 


0.92 


0.95 


0.75 


0.85 


0.95 


0.70 


0.47 


0.84 


0.18 


1.14 


2.50 


1.15 


0.78 
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Figure 1: Time dependence of A^, Xp, Xk and Xk* without the source terms (sohd curves), 
with the quark-interchange, annihilation and resonant processes (dashed curves) and with 
only the quark- interchange processes (dotted curves). 
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